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Thermodynamic and transport properties of normal disordered conductors are strongly influenced 
by the proximity of a superconductor. A cooperation between mesoscopic coherence and Andreev 
scattering of particles from the superconductor generates new types of interference phenomena. We 
introduce a field theoretic approach capable of exploring both averaged properties and mesoscopic 
fluctuations of superconductor/normal- metal systems. As an example the method is applied to the 
study of the level statistics of a SNS-junction. 



Physical properties of both superconductors and meso- 
scopic normal metals are governed by mechanisms of 
macroscopic quantum coherence. Their interplay in hy- 
brid systems comprised of a superconductor adjacent to a 
normal metal gives rise to qualitatively new phenomena 
(see Ref. § for a review) : Aspects of the superconducting 
characteristics are imparted on the behaviour of electrons 
in the normal region. This phenomenon, known as the 
"proximity effect" , manifests itself in a) the mean (dis- 
order averaged) properties of SN-systems being substan- 
tially different from those of normal metals and b) vari- 
ous types of mesoscopic fluctuations which not only tend 
to be larger than in the pure N-case but also can be of 
qualitatively different physical origin. Although powerful 
quasi-classical methods, based largely on the pioneering 
work of Eilenberger Q and Usadel (^] , have been devel- 
oped to analyse the manifestations of the proximity effect 
in average characteristics of SN-systems, far less is known 
about the physics of mesoscopic fluctuations: While the 
quasi-classical approach is not tailored to an analysis of 
fluctuations, standard diagrammatic techniques used 
in the study of N-mesoscopic fluctuations can often not 
be applied due to the essentially non-perturbative influ- 
ence of the fully established proximity effect. Important 
progress was made recently by extending the scattering 
formulation of transport in N-mesoscopic systems to the 
SN-case Q. This approach has proven powerful in the 
study of various transport fluctuation phenomena but is 
not applicable to the study of fluctuations on a local and 
truly microscopic level. 

In the present Letter we introduce a general frame- 
work that combines key elements of the quasi-classical 
approach with more recent methods developed in N- 
mesoscopic physics into a unified approach. As a result 
we obtain a formalism that can be applied to the gen- 
eral analysis of mesoscopic fluctuations superimposed on 
a proximity effect influenced mean background. In order 
to demonstrate the practical use of the approach we will 
consider the example of spectral fluctuations as a typical 
representative of a mesoscopic phenomenon. The density 
of states (DoS) of N-mesoscopic systems exhibits quan- 
tum fluctuations around its disorder averaged mean value 
which can be described in terms of various types of uni- 



versal statistics. The analogous question in the SN-case - 
What types of statistics govern the disorder induced fluc- 
tuation behaviour of the proximity effect influenced DoS? 
- has not been answered so far. Although space limita- 
tions prevent us from discussing this problem in detail, 
our main result, the emergence of some kind of modified 
Wigner Dyson statistics ||, will be derived below. 




FIG. 1. Geometry of the SNS-junction. 

To be specific we consider the geometry shown in 
Fig. [l], where tp — (pi — tp2 represents the relative phase 
of the complex order parameters of the two superconduc- 
tors. It is well known Q that even the mean DoS of the 
SNS-system exhibits non-trivial behavior which is diffi- 
cult to describe within standard perturbative schemes: 
States which fall within the superconducting gap, A, 
are confined to the normal metal. The proximity ef- 
fect then further induces a minigap in the DoS of the 
normal region around the Fermi-energy, tp, whose size 
of 0(E C = D n /L 2 ) depends sensitively on ip (D n is the 
diffusion constant and fi — 1 throughout). To analyse 
the fluctuation behavior of the DoS, v(e), around its dis- 
order averaged background, {v(£j), we will consider the 
two-point correlation function 
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The starting point of our analysis is the Gor'kov equa- 
tion for the matrix advanced/retarded (a/r) Green func- 
tion fll 
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c± = e ± iO, A is the vector potential of an external mag- 
netic field, A = Acr ph exp(— iipa^ 1 ) represents the (spa- 
tially dependent) complex order parameter with phase tp, 
and Pauli matrices cr ph operate in the Nambu or parti- 
cle/hole (ph) space. The impurity potential in the N- 
region is taken to be Gaussian (^-correlated with zero 
mean and correlation {V(r)V{f)) = 5 d (r — r 1 ) /(2ttvt), 
where v denotes the DoS of the bulk normal metal at ep, 
and t represents the mean free scattering time. In the 
following the complex order parameter in the S-region 
is imposed and not obtained self-consistently |^]. Where 
the S- and N- region are distinct (as in the SNS junc- 
tion), the bulk DoS, v n>s and scattering time, T n>s will be 
chosen independently. 

Traditionally the impurity averaged Green function (|^) 
is computed within a quasi-classical approximation, i.e. 
the Schrodinger equation (|^) is reduced to an effective 
transport equation, the Eilcnberger equation [g], which 
in the dirty limit simplifies further to the diffusive Us- 
adel equation || . Here we develop a field theoretical for- 
mulation that integrates concepts of the quasi-classical 
formalism into a more general framework allowing for 
the computation of disorder averaged products of Green 
functions, a necessary requirement for the calculation of 
correlation functions such as (Q). The basic strategy will 
be to start from a (microscopically derived) generating 
functional whose points of stationary phase obey the Us- 
adel equation. By investigating fluctuations around this 
quasi-classical limit, correlations between the different 
Green functions will be explored. In the following we 
formulate this program in more detail. 

As in the pure N-case, ensemble averaged products of 
advanced and retarded Gorkov Green functions can be 
described in terms of generating functionals of nonlinear 
cr-model type |n| (see Ref. |j) for a review on the cr-model 
analysis of Green functions in N-mesoscopic physics). In 
the dirty limit, (e, A) < r _1 <C ep, the generalization of 
the supersymmetric N-type cr-model M reads 
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where d = d — i(e/c)[A^al r <g> cr ph , -J represents a covari- 
ant derivative, A<p = A + cj (2e)d(j> accounts for both 
the external field and the phase of the order parameter, 
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the Pauli matrices ct"\ cr tr and cr a 



oper- 



ate in fermion/boson, time-reversal and ar-blocks respec- 
tively ||. The symbol D stands for a space dependent 
diffusion constant which may take separate values, de- 
noted as D nyS , in the N and S regions. Although spe- 
cific pre-exponential source terms (denoted by ellipses 



in Eq. (Q)) must be chosen according to any given cor- 
relation function (such as (0)), their precise form does 
not influence the analysis below and we therefore refer 
to Ref. H for their detailed structure. The integration 
in (|J) extends over a 16 x 16-dimensional matrix field 
Q = T _1 CT3 h i8)cr| r T, whose symmetries are identical with 
those of the conventional cr-model || . 

The expression (|J) differs in two respects from the cr- 
model for N-systems: i) the appearance of a ph-space 
associated with the 2 x 2-matrix structure of the Gorkov 
Green function, and ii) the presence of the order param- 
eter A. Whereas i) can be accounted for by a doubling 
of the matrix dimension of the field Q, ii) calls for more 
substantial modifications: For A =/= standard pertur- 
bative schemes for the evaluation of the functional (Q) 
fail |ll| , an indication of the fact that the superconduc- 
tor influences the properties of the normal metal heavily. 
Under these conditions a more efficient approach is first 
to subject the action to a mean field analysis and then 
to consider fluctuations around a newly defined — and 
generally space dependent — stationary field configura- 
tion. A variation of the action (||) with respect to Q, 
subject to the constraint Q 2 — 1, generates a non-linear 
equation for the saddle-point, 



Ddi(QdiQ) + |Q, Aaf ~ i(e + w + CT 3 ar )CT 3 pn j = 0. (5) 

Current conservation implies the boundary condi- 
tion Q, a n Qd x \ Xn Q = a s Qd x \ Xa Q, where ct„ jS = 
e 2 v n ,sD n . s denotes the conductivity and d x \ Xn , . is a nor- 
mal derivative at the N(S)-side of the interface. 

An inspection of (||) shows that only the particle/hole 
components of the matrix field Q are coupled by the sad- 
dle point equation. It is thus sensible to make a block 
diagonal ansatz Q = bdiag (q+ , g_) , where the eight di- 
mensional retarded, q+ and advanced, g_ subblocks are 
diagonal in both time reversal and boson/fermion space. 
Noting that the saddle point configuration —i-nvq± of the 
nonlinear cr-model is associated with the impurity aver- 
aged retarded/advanced Green function B, we identify 
Eq. (||) as the Usadel equation. The general connection 
between the cr-model formalism and the quasi-classical 
approach has first been noticed in p2| . 

Eq. (^|), in its interpretation as the Usadel equation, 
has been discussed at length in the literature |^, , ff5[ . 
Although complex in general, the solutions have a sim- 
ple qualitative geometric interpretation: Employing the 
explicit parametrization q± = q± ■ cr ph , Eq. (||) describes 
the gradual rotation of the three dimensional vector q 
from a direction almost parallel to e\ in the bulk super- 
conductor to a value aligned with £3 deep in the normal 
metal. Closed analytical solutions have been obtained for 
various types of geometries, including the SNS-structure 
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above. 

So far our analysis has been for SN-systems of a general 
geometry. Specializing the discussion to the particular 
SNS-junction shown in Fig. (|l|), we set A(r) = A6(|x| — 
L/2) constant inside the superconductor (A 3> E c ), and 
zero in the normal region, with a phase n/2 + sgn(x)p>/2. 
The saddle-point equation depends sensitively on both 
the presence or absence of an external magnetic field and 
the phase difference between the order parameters. Tak- 
ing the external field to be zero, it is convenient to fo- 
cus on two extreme cases: (i) ip = (orthogonal sym- 
metry), and (ii) p ^> l/ifg (unitary symmetry). Here 
g = E c /d ^S> 1 denotes the dimensionless conductance 
and d represents the bulk single-particle level spacing of 
the normal metal. 

The disorder averaged local DoS can be obtained from 
the analytical solution of the Usadel-saddle point equa- 
tion |^,[I^-|l7| as v(r) — vRe{q + (r)] 3 . Its space/energy 
dependence is shown in Fig. ||for the case (i), zero phase 
difference, and a particular value of the material param- 
eter 7 = v n \JD n jv s \JD s . The most striking feature of 
the average DoS is the appearance of a spatially constant 
minigap in the N-region. The gap attains its maximum 
width E c at ip = and shrinks to as p approaches 
ir [^5|. Within the superconductor, v{E c < e <C A) de- 
cays exponentially on a scale set by the bulk coherence 
length £ = (D s /2A) 1 /2. 




FIG. 2. Local DoS v shown as a function of energy and 
position for 7 = 0.1, L/t; — 5 and ip — 0. 

We next turn to the main subject of this Letter, the 
issue of fluctuations around the Usadel saddle-point. Em- 
ploying the parametrization Q = T~ 1 QT 1 T ^ 1, three 
qualitatively different types of fluctuations can be iden- 
tified: (a) fields T which are diagonal in the space of 
advanced and retarded components, (b) T's which com- 
mute with all matrices af but mix advanced and re- 



tarded components, and (c) T's fulfilling neither of the 
conditions (a) and (b). Fluctuations of (a)-type preserve 
the ar-diagonal structure of the saddle point. These fluc- 
tuations do not give rise to correlations between advanced 
and retarded Green functions. Nonetheless, they are of 
physical significance: Quantum corrections to the Us- 
adel solution, most importantly the renormalization of 
the minigap by weak localization effects and the exis- 
tence of rare prelocalized states ]i~^ , |l8| below the gap, 
are described by fluctuations of this type. For brevity, 
we postpone further discussion of these results to a sep- 
arate paper JfTj and, instead, turn to the discussion of 
the second type of fluctuations, (b). 

Consider the saddle point equation (^]) in the simple 
case u) = p = 0. Obviously, as it commutes with all ma- 
trices cr? h , any spatially constant rotation T of type (b) 
gives rise to another solution. In other words, the (b)- 
fluctuations represent Goldstone modes with an action 
that vanishes in the limit of spatial constancy and uj — > 0. 
Since any T diagonal in ph-space inevitably has to couple 
between advanced and retarded indices ||, these modes 
lead to correlations between advanced and retarded Green 
functions (and thereby to mesoscopic fluctuations) which 
become progressively more pronounced as ui approaches 
zero. 

In the limit of small frequencies u> < E c , the ergodic 
regime, the global zero mode Qo = T _1 gTo, [To,a ph ] = 
0, To(r) — const., plays a unique role: Whereas fluctua- 
tions with non-vanishing spatial dependence give rise to 
contributions to the action of 0(g ^> 1) ||, this mode 
couples only to the frequency difference uj. Restricting 
attention to the pure zero mode contribution, we obtain 
the effective action 

SolQo] =-i~^stT[Q a?], (6) 
2 d(e) 

where d(e) = (J v(e))~ 1 denotes the average level spac- 
ing and the ph-degrees of freedom have been traced out. 
From this result it follows H that, in the ergodic regime, 
the spectral statistics of an SNS system is governed by 
Wigner-Dyson fluctuations (|Jl^l of (i) orthogonal or (ii) 
unitary symmetry superimposed upon an energetically 
non-uniform mean DoS. Furthermore, a comparison of 
Eq. (H) with the analogous action for N-systems H shows 
the correlations to depend on an average level spacing 
that is effectively halved. This reflects the strong "hy- 
bridization" of levels at energies ~ epit induced by 
Andreev scattering at the SN-interface. 

In further contrast to N-systems, the range over which 
Wigner-Dyson statistics apply turns out to be greatly 
diminished by non-universal fluctuations. To prove the 
last statement, fluctuations of type (c), coupling between 
advanced/retarded and particle/hole components simul- 
taneously, have to be taken into account. Due to their 
complexity, the detailed analysis of the (c)-type fluctua- 
tions is cumbersome and will be deferred to a forthcom- 
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ing publication |l7| . Here we restrict ourselves to a brief 
and qualitative discussion of the principal effect of these 
fluctuations on the universal zero mode action @. 

As can be seen from the general structure of the ac- 
tion (0), (c)-type fluctuations in the vicinity of the mini- 
gap are generally 'massive' (governed by an action which 
is at least of order e/d > g ^> 1). It is thus permissible to 
treat these fluctuations in a simple Gaussian approxima- 
tion. As a result we arrive at an effective action which 
is non-perturbative in the above zero mode configuration 
and quadratic in the (c)-type perturbations. Integrating 
over these fields in a spirit similar to the analysis per- 
formed in Ref. p(if we obtain the modified zero mode 
action 
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where fc(e) ~ O(l) denotes a constant dependent on the 
sample geometry p7| . Eq. (|t]) has a structure equiva- 
lent to that found in the study of universal parametric 
correlation functions and explicit expressions for i?2 for 
both orthogonal and unitary ensembles can be deduced 
from Ref. [|TJ . Qualitatively, the additional contribution 
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vanishing frequencies to. Already for energy separations 
u>/d(e) ~ the zero- mode integration is largely sup- 
pressed which manifests in an exponential vanishing of 
the level correlations on these scales. This is in contrast 
to the pure N-case where the Wigner- Dyson regime (pre- 
vailing up to frequencies u> ~ E c ) is succeeded by other 
forms of algebraically decaying spectral statistics in the 
high frequency domain ui > E c |22| . 

In conclusion a general framework has been developed 
in which the interplay of mesoscopic quantum coherence 
phenomena and the proximity effect can be explored. An 
investigation of the spectral statistics of an SNS geom- 
etry revealed that level correlations are Wigner-Dyson 
distributed with strong non-universal corrections at large 
energy scales. Finally, we remark that for quantum 
structures in which transport is not diffusive but bal- 
listic and boundary scattering is irregular, a ballistic o- 
model involving the classical Poisson bracket can be de- 
rived |23). In this case, the saddle-point condition recov- 
ers the Eilenberger equation of transport j!| . 
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